function pressure_dH = fn_int_pressure_dH(vx, vep, vep_dH, T)

vx_nt = fn_nt(vx);
pressure_dH = -2*sum(vx_nt.*vep_dH.*fn_G(vep, T))/2/pi;

end